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ABSTRACT 


The  various  design  considerations  of  a  solar  ephaneris  synthesizer  for  the 
LES-8/9  stationkeeping  system  are  discussed  here.  Different  approximations  are 
used  to  reduce  the  amount  of  hardware  required  for  building  the  synthesizer. 

The  overall  accuracy  of  the  synthesizer  is  about  0.01  deg  for  orbits  inclined 
1:5)  to  5.0  deg  from  the  ecliptic  plane.  For  orbits  inclined  at  greater  angles 
iq)  to  10  deg,  the  short  term  accuracy  is  0.01  deg/day  and  the  long  term 
accuracy  is  about  0.05  deg/year.  A  few  ephemerides  generated  by  the  synthe¬ 
sizer  are  graphically  illustrated. 
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DESIGN  OF  A  SOLAR  EPHEMERIS  SYNTHESIZER  FOR  LES-8/9 
I.  INTRODUCTICN 

An  autonomous  stationkeeping  system  for  a  synchronous  satellite  is 
required  to  control  the  mean  longitude  of  the  satellite.  The  daily  variation 
of  the  longitude  about  a  mean,  which  results  from  the  eccentricity  and 
inclination  of  the  orbit  to  the  equator,  is  uncontrolled  and  goes  through  a 
well  predictable  cycle.  The  input  required  for  the  above  stationkeeping 
‘  ystem  is  the  angular  position  of  the  satellite  in  its  orbit  plane.  When 
the  eccentricity  of  the  orbit  is  not  required  to  be  controlled,  only  two 
measurements  at  opposite  points  in  the  orbit  are  sufficient.  For  a  synchro¬ 
nous  satellite  in  circular  orbit,  the  anpalar  position  can  be  determined  by 
observing  the  azimuthal  conponent  of  the  sun’s  motion  from  a  perfectly  oriented 
satellite,  that  is,  a  satellite  having  no  attitude  errors.  The  sun’s  motion 
need  not  be  observed  continuously,  but  it  is  sufficient  to  know  the  time  of 
sun  transit  across  two  azimuth  transit  sensors  looking  in  opposite  directions. 
Time  is  kept  accurately  by  a  stabilized  frequency  source  having  a  short  term 
stability  of  0.5  sec/day  and  a  long  terra  stability  of  2.5  sec/year.  Thus, 
the  orbit  angle  of  the  satellite  can  be  computed  at  the  time  of  sun  transit 
by  comparing  the  observed  transit  time  with  the  expected  transit  time.  The 
expected  transit  time  is  a  function  of  inclination  and  ascending  node 
of  the  orbit.  The  mean  interval  between  the  transit  times  of  the  sun  across 
a  transit  sensor  over  one  year  is  24  hours.  The  difference  between  the 
mean  transit  time  and  expected  transit  time  is  a  function  having  a 
period  of  one  year.  This  function  can  alternatively  be  ejqiressed  as  an  angle 
rather  than  time  with  the  usual  conversion  factor  of  360°/24  hr.  The  solar 
ephemerides  for  observation  from  a  perfectly  oriented  satellite  are  usually 
expressed  as  angular  deviations  at  regular  time  intervals,  vhereas  the  above 
function  gives  the  time  difference  between  the  expected  and  mean  transit  times 
at  the  expected  transit  times.  In  other  words,  the  ephemerides  are  samples  of 


angular  deviations  at  uniform  time  intervals,  whereas  the  function  considered 
here  is  a  sanple  of  the  time  delay  at  non-uniform  time  intervals.  The  function 
when  expressed  in  angular  units  will  be  indistinguishable  from  the  solar  ephem- 
eris  within  our  desired  accuracy.  Therefore,  without  making  any  further  dis¬ 
tinction  we  shall  refer  to  the  time  function  expressed  ir^  angular  units  as  the 
solar  ephemeris  and  a  synthesizer  for  the  function  as  an  ephemeris  synthesizer. 

In  the  above  discussion  it  was  assumed  that  the  satellite  is  perfectly 
oriented.  This  is  true  only  within  a  fraction  of  a  degree  for  a  three  axis 
attitude  controlled  satellite  similar  to  LES-8/9.  Therefore,  the  satellite 
serves  as  a  non-staticMiary  platform.  However,  the  correction  for  suri  observa¬ 
tion  can  be  conveniently  ^lied  for  small  deviations  of  the  satellite  atti¬ 
tude  [1] .  The  error  due  to  eccentricity  can  be  averaged  out  to  a  high  degree 
by  taking  the  average  of  an  even  number  of  consecutive  observations  taken  at 
opposite  ends  of  the  orbit. 

For  the  LES-8/9  satellites,  viiich  will  be  nominally  in  the  ecliptic  plane, 
it  is  sufficient  to  be  able  to  generate  ephemeris  for  orbits  inclined  up  to 
ten  degrees  from  the  ecliptic  plane.  The  short  term  accuracy  of  the  synthes¬ 
ize;  is  desired  to  be  good  enough  to  give  a  high  fuel  efficiency  for  the 
stationkeeping  system.  For  a  ±0.1  deg  stationkeeping  system  a  desirable 
value  of  the  short-term  accuracy  is  about  0.01  deg  [2].  The  long  term 
error,  which  is  not  veiy  much  affected  by  the  averaging  scheme  and  has  a  very 
small  effect  on  the  fuel  efficiency  catises  a  slow  drift  of  the  satellite.  For 
orbits  having  inclination  up  to  5.0  deg  from  the  ecliptic,  the  desired  long 
term  accuracy  is  0.01  deg.  For  orbits  at  greater  inclinations  the  long  term 
error  could  be  allowed  to  be  larger  at  the  cost  of  a  slow  station  drift.  The 
limiting  case  of  10  deg  inclination  the  drift  is  about  0.05  deg/year. 

This  tolerance  allows  a  sinple  design  and  a  small  overall  size  for  the 
synthesizer. 

In  this  note  the  various  design  considerations  for  the  synthesizer  are 
discussed.  Solar  ephemerides  synthesized  for  a  few  different  orbits  by  the 
synthesizer  are  illustrated  in  the  appendix. 
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II.  NOTATIONS 


The  following  notations  are  used  in  this  note, 
r  =  radius  vector  from  the  satellite  to  the  sun 
y  =  GM  ,  where  G  is  the  gravitational  constant  and  M  is  the 
mass  of  the  earth 
a  =  semi -major  axis  of  the  orbit 
e  =  eccentricity  of  the  orbit 

f  =  true  anomaly  of  the  orbit 

n  =  mean  angular  velocity 

i  =  mean  anomaly 

u  =  argument  of  the  perigee 
=  Bessel  function  of  order  n 
U  =  0)  +  f 


a,, a 


1»“2*** 


t 


x,y,2 


Q 

i 

0 


eccentricity  coefficients  in  the  expansion  of  the  true  anomaly 
as  a  function  of  mean  anomaly 
time  in  days 

the  difference  between  true  anomaly  and  mean  anranaly  for  the 

ecliptic  plane 

transpose  of  a  vector 

rectangular  coordinates  of  a  vector 

right  ascension  of  ascending  node 

inclination  of  the  orbit 

error  in  position  measurement  after  averaging  measurements  at 
two  opposite  ends  of  the  orbit 
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inclination  coefficient 


.2  .4 

-  ^  ^  - 

il»2Cu)  =  part  of  a  solar  ephemeris  resulting  from  the  inclination  of  the 
satellite  orbit 

ijjjCu)  =  an  error  in  ^roximation 

=  an  error  due  to  solar  parallax  from  synchronous  orbit 
ipgCk)  =  an  error  due  to  junp  summation. 

Int(s)  =  integral  part  of  "s" 

a  -  angle  of  the  sun-earth  vector 
6  =  earth  satellite  angle 
T  =  2-rr  27  ^  ,  a  time  variable 
Y  =  earth-sun  angle  subtended  at  thr  satellite 
Ti  =  inclination  of  the  equatorial  plane  with  the  ecliptic  =  23.44° 
e^  =  eccentricity  of  satellite  orbit 
*^ij  ~  Kroneker  delta,  1  for  i  =  j  ,  0  for  i  5^  j 


III.  SUN  AZIMJTH  FOR  ECLIPTIC  ORBIT 

For  our  desi’.ed  accuracy,  the  motion  o£  a  satellite  around  the  sun  may  be 
approximated  by  the  motion  of  the  earth  around  the  sun.  The  small  inaccuracy 
in  this  approximation  may  be  treated  as  a  parallax  error.  It  is  assumed  that 
the  satellite  is  in  a  circular  orbit  and  perfectly  oriented,  so  that  its  equa¬ 
torial  plane  is  coplanar  with  the  orbit  plane. 

The  motion  of  the  earth  arotuid  the  sun  may  be  treated  as  a  two  body 
classical  problem.  Therefore, 

yv 

$  r  f., 

r  =  -  y  ,  Ci- 

T 

where  r  is  the  sun-earth  vector  and  p  is  the  product  of  the  gravitational 
constant  and  the  mass  of  the  sun. 

The  solution  of  the  above  equation  is  the  ellipse  [3]* 

r  =  „  (3- 

1+e  cos  £  *  ^ 

where  a  is  the  semi -major  axis,  e  is  eccentricity  and  f  is  true  anomaly 
of  the  earth  orbit. 

The  mean  angular  velocity  n  is  given  by 

1/2  -3/2 

n  =  y  '  a  .  (3- 

The  mean  anomaly  Si  is  given  by 

i  =  nt  ,  (3- 


where  t  is  the  time  measured  from  the  initial  perigee  passage.  A  series 
solution  of  Eq.  (3-1)  is  given  by  [3]. 


£  =  £ 


00  00 

*  ■*  I 

S=1  '  P=1  ^  ' 


The  number  vrithin  the  square  bracket  []  refers  to  the  reference  nmiber. 
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where  J  is  Bessel's  fianction  of  order  s,  and 
s 

6  =  .  (3-6) 

Equation  (3-5)  written  more  exolicitly  gives  [3] 

f  =  £  +  ^2e  —  ^  e^  +  ^  ^  *  *  *)  ^ 

+  .  .j  sin  4£  +  e^  -  . .  .^  sin  S£ 

^  •  •  •  « 

=  £  +  sin  £  +  a2  sin  2£  +  a^  sin  3£  +  a^  sin  4£  +  ...  ,  (3-7) 

where  a2,  a^,  ...  are  the  coefficients  of  the  eccentricity  terms  and 

271 1 

£  =  27  •  Therefore  the  variation  of  the  sun  azimuth  angle  from  the  mean 

sun  angle  is  given  by 

ip^(£)  =  f  -  £ 

=  sin  £  +  32  sin  2£  +  a^  sin  3£  +  a^  sin  4£  +  ...  (3-8) 

Notice  that  tj^-,(£)  is  an  odd  function  around  £  =  nv,  n  =  o,±l,±2,... 
but  it  is  not  an  even  function  around  £=(n  +  j)TT,  n  =  o,  ±1,  ±2,.. 

therefore,  the  sign  of  the  function  cannot  be  changed  by  changing  the  phase  by  a 
an  angle  ir  . 


IV.  SUN  AZIMJTH  FOR  NON-EaiPTIC  ORBITS 

The  satellite's  plane  of  motion  is  assumed  to  be  the  XT  plane  in  Fig.  1. 
OS  is  the  satellite-sun  vector.  The  sun  azimuth  is  obtained  by  projecting 
the  vector  OS  on  the  XY  plane.  The  coordinates  of  the  sun  in  XYZ  frames 
are  obtained  bv  the  following  transformation 


A 
s 

I  ^ 


(1) 

Rotating 

[r,0,0]^ 

by  -  u 

around 

1" 

(2) 

Rotating 

[r’, 

O 

o 

by  —  i 

around 

X' 

(3) 

Rotating 

[r”, 

o,of 

by  -  12 

around 

z 

pos  fl 

sin 

o' 

1 

0 

0 

cos  u 

sin  u 

-sin  cos  12 

0 

0 

cos  i 

sin  i 

-sin  u 

cos  u 

L» 

0 

1 

0 

-sin  i 

cos  i 

0 

0 

r” 

0 

0 

J 

m  « 

r  cos  u  cos  —  r  sin  u  sin  n  cos  i 
r  cos  u  sin  fl  +  r  sin  u  cos  cos  i 
r  sin  u  sin  i 


(4-1) 


For  convenience  let  the  X  axis  be  coincident  with  the  ascending  node, 
i.e.,  =  0  ,  Eq.  (4-1)  becomes: 


;  1 

-  - 

X 

r  cos  u 

♦  -aw 

i  1 

Y 

= 

r  sin  u  cos  i 

1 

(4-2)  i  1 

«  ^ 

Z 

r  sin  u  sin  i 

;  1 

1  % 

i  i 


I  ^ 

I  ai 

i  I 
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From  Eq.  (4-2),  the  sun  azimuth  is  given  by 


0 


.  -1  sin  u  cos 
=  sin 


PUJ 


(4-3) 


where 

...  L  -  2  .2  .M/2 

p(i)  =  f  1  —  sin  u  sin  ij  ' 

The  difference  between  the  sun  aziimith  angles  of  an  ecliptic  orbit  and  a 
non-ecliptic  orbit  is  given  by 


ij/2(u)  =  u  -  9 

(4-4) 


(4-5) 

From  Eq.  (4-4),  i|»2(u)  can  be  shown  to  be  extremum  at 


_  — sin  2u/ 

'  -pnri 

1.,  ..  sin  2u  .  1  .. 

=  j(l  -  cos  1)  *  ^(1  -  cos  I) 

=  A-  ^  sin  2u  —  ^  sin  4u  +  . . . 


3 

..3  sin  2u 
P(i)^ 
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V.  i\PPROXIMA.TION  OF  :HE  SUN  AZIMjni 

The  inclination  correction  term  for  the  orbit  given  by  Eq.  (4-5)  is  not 
ver>'  convenient  for  building  a  siii5?le  solar  ephemeris  synthesizer.  Let  i{'2(u) 
be  approximated  by  the  first  term  of  Eq.  (4-5).  Therefore,  the  approximate 
!|>2(u)  is  given  by 

'  fr  *  y 


=  I(i)  sin  2u  , 


(5-1) 


where 


The  values  of  I (i)  are  tabulated  for  various  inclir.itions  in  Table  T . 

_ 

The  maxiimim  error  in  ij/'2(u)  approximation  is  about  ,  and  it  is  given 
in  Table  II  for  various  inclinations. 

From  Eq.  (3-7) 


I(i)  sin  9,  =  j^j^(^)  —  32  sin  ZZ  — 


sin  3Jl  —  a^  sin  4£.  . . .  .j 
(5-2) 


Let  £  =  2u  ,  then 


^  [ij(2u)  - 


a-  sin  4u  —  a,  sin  6u  —  a.  sin  8u  + _ 

2  3  4 


=  llil 


'j'j^(2u)  —  1^2  sin  4u  +  a^  sin  6u  +  a^  sin  8u  +  ...^ 


(5-3) 


TABLE  I 

TABLE  OF  THE  INCLINATION  FACTOR  I(i) 


Inclination  i 

deg 

I(i) 

deg 

0 

0 

2.5 

0.0271 

5 

0.109 

10 

0.435 

20 

1.739 

23.5 

2.33 

TABLE  II 

TABLE  OF  THE  ERROR  IN  ^(u)  APPROXIMATION 

DUE  TO  NEGLECTING  FOURTH  AND  HIGHER  HARMONICS 

Inclination  "i" 

Error  = 

deg 

deg 

0 

0 

2.5 

0.0004 

5 

0.0016 

10 

0.0066 

20 

0.0264 

23.5 

0.0501 

11 


Let  approximated  by  the  first  term  in  Eq.  (4-8).  The  maximum  error 

’^,(u)  in  this  approximation  satisfies 

’^3^^-’  ^  ^a^  1^3!  I34I  "^  ••••]  • 

For  the  earth  orbit  e  =  0.01675,  and  the  eccentricity  terms  in  Eq.  (3-7)  are 

a^  =  3.350  X  10"^  , 

a2  =  3.506  x  lo""^  , 

a^  =  5.09  X  10"^  ,  ...  .  (5-5) 

Using  the  above  in  Eq.  (5-4),  we  get 

V  .(u)  £  I(i)  X  1.05  X  10’^  rad 

<  0.603  I(i)  deg  .  (5-6) 


Table  Til  gives  the  maximum  value  of  tf»2(u)  for  various  inclinations.  The 
overall  variation  of  the  sun  azimuth  due  to  the  earth's  eccentric  orbit  and 
the  inclination  of  the  satellite's  orbit  plane  with  the  equator  is  approximated 
by 

’^(t)  = 


(5-7) 


where  u  is  the  argument  of  the  perigee  (see  Fig.  1). 

In  the  above  discussion  the  effect  of  parallax  due  to  observations 
from  the  synchronous  orbit  has  not  been  considered.  The  small  parallax  cor¬ 
rection  is  necessary  because  in  the  Eq.  (5-7)  the  synchronous  distance  is 
assumed  to  be  negligible  compared  to  the  sun-earth  distance.  The  maximum 
solar  parallax  observed  from  a  synchronous  orbit  in  the  ecliptic  plane  is: 
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TABLE  III 


TABLE  OF  MAXIMUM  ERROR  tp^Cu)  IN  THE  ijijCu)  APPROXIMATION 


i 

deg 

deg 

0 

0 

2.5 

0.0003 

5 

0.0011 

10 

0.0046 

20 

0.0182 

23.5 

0.0252 

TABLE  IV 

maxim:.!  error  in  inclination  correction  term 

DUE  TO  JUMP  SUNWATION 


i 

deg 

a^/I 

’^5 

deg 

0 

00 

0 

2.5 

70.4 

0.00002 

5.0 

17.6 

0-003 

6.6 

10.0 

0.010 

10.0 

4.4 

0.050 

14.8 

2 

0.240 

r,  _  ^^syn  _  2  X  4.2165  x  lo"^ 

-  A.U.  '  1^496,1.8 

=  0.03230  deg  ,  (5-8) 

where  is  the  synchronous  radius  and  A.U.  is  the  actronondcal  unit 

(mean  earth-sun  distance).  The  solar  parallax  is  almost  constant,  if  the 
measurement  for  the  sun  azimuth  is  performed  for  the  same  earth-sun  angle 
subtended  at  the  satellite.  This  constant  error  can  be  eliminated  by  either 
taking  the  average  of  measurements  at  two  opposite  points  in  the  orbit  or 
by  adjusting  the  alignment  of  sun  sensors. 

Another  source  of  error  in  using  Eq.  (5-7)  for  the  synthesis  of  solar 
ephemeris  is  due  to  discrete  storage  of  the  function  4;^(t).  Differentiating 
Eq.  (3-8) 

%1\27  pi  ^  *  ^®^3  •  *  •] 

>  ^2^  ^3*  •••  from  (5—5)  xn  Eq.  (5—9) 

'*'l^*'^]max  "  365.27  Pl  ^^2  ^3  “) 

=  0.0338  6t  deg  .  (5-10) 

If  the  function  ij»j^(£)  be  sanpled  at  j  day  intervals  then  the  maximum 
difference  between  the  saaiqjling  intervals  is  0.017  deg,  and  the  maximum 
error  due  to  discrete  storage  is  0.0084  deg. 


6  4'i(^) 
~Zt 

Using  the  constants 
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VI .  l-MORS  DUn  TO  ORBIT  PIIRTURBATIONS 

The  orbit  of  a  synchronous  satellite  changes  significantly  over  its  life¬ 
time  due  to  geopotential,  luni-solar  potential  and  solar  radiation  pressure 
[4,5,6].  The  effect  of  other  perturbations  are  small.  The  solar  ephemeris  for 
the  satellite  is  affected  by  the  changes  in  orbital  elements.  As  an  example 
consider  a  synchronous  satellite  in  an  eccentric  orbit  in  the  ecliptic  plane 
(Fig.  2) .  From  liq.  (3-7)  the  angle  of  the  sun-earth  vector  01:  is  given  by 

or) 

a  =  SL  —  ^  ^  aj^  sin  k(2.  —  £.^)  ,  (6-lJ 

k=l 

where  is  the  mean  anomaly  for  the  direction  OP  and  £  =  3^5-27  * 
arly,  the  earth  satellite  angle  SE  is  given  by 

00 

6  =  T  -  Tq  sin  j(T  -  t^)  ,  (6-2) 

where  is  the  angle  between  the  earth  earth-satellite  vector  when  the  satel¬ 
lite  is  at  its  perigee  and  the  initial  sun-satellite  vector  OP.  The  variable 

77t' 

T  =  2Tr  and  b^ 's  are  the  coefficients  similar  to  a^'s  when  the 

eccentricity  e  is  substituted  by  the  eccentricity  of  the  satellite  orbit  e^  . 
Frcmi  (6-1)  and  (6-2),  the  earth-sun  angle  at  the  satellite  is  given  by 


Y  =  6  —  a 


=  (T  -  1)  *  fT„  - 

k=l 


aj^  sin  k(Il  -  y  +  bj^  sin  k(T  -  t^) 


^"^o  ^o^  *  365.27 


CO 

E 

k=l 


2Trk 


366.27 


365.27  V  \  365.27 


(6-3) 
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If  it  is  assumed  that  the  satellite  has  no  pi..ch,  roll  or  yaw  error,  then 
the  successive  times  of  the  sun  coincidence  by  a  sun  azimuth  sensor  on  the 
satellite  are  obtained  by  solving  (6-3)  for  Y  =  ^  2mr,  n-  0:1,2,..., 

is  the  angle  between  the  azimuth  sensor  field  of  view  and  the  nominal  pitch 
direction  of  the  satellite.  A  general  series  solution  or  a  closed  form  solu¬ 
tion  of  Eq.  (6-3)  is  a  non- trivial  problem.  However,  it  is  seen 
in  (5-5)  tiiat  aj^'s  converge  to  zero  very  rapidly.  Similarly  even  for  an 
eccentricity  of  0.01  the  bj^'s  will  converge  to  zero  quite  fast.  To 
eliminate  the  effect  of  the  variation  due  to  the  eccentricity  of  the  satellite's 
orbit,  the  sun  azimuth  angle  is  measured  at  opposite  ends  of  the  orbit,  that 
is  at  Y  =  Yq  +  Znir  and  Y  =  Yq  +  (2n  +  1)  it  ,  n  =  o,  1,  2,  ...  The  maxiraun 
error  in  the  position  ruiasurement  of  the  satellite  when  the  average  of  two  such 
measurements  are  tai;en  is  given  by 


0 


max 


+  -hJs 

365.27  365.27/ 


4ire„ 

5^  _ 

365.27  ’ 


(6-4) 


where  e^  is  the  eccentricity  of  the  satellite  orbit,  ro'  an  orbit  having 
Sg  =  0.01  ,  this  error  amounts  to  0.01  deg.  This  error  is  of  a  long  term 
period  of  one  year  and  therefore  quite  insignificant.  For  an  inclined  orbit 
the  averaging  scheme  will  eliminate  the  eccentricit>’  error  in  the  sane  manner. 
The  residual  error  after  averaging  is  the  sane  order  as  (6-4)  for  orbits  having 
inclination  vp  to  10  deg.  The  position  measurement  scheme  with  averaging  is 
therefore  not  sensitive  to  the  variations  of  eccentricity  and  argument  of  the 
perigee. 

The  effect  of  the  changes  in  the  inclination  and  the  ascending  node  can  be 
conputed  from  Eq.  (5-7).  The  variation  due  to  changes  in  these  two  elements  is 

M(t)  =  61  .  |S  '  (6-5) 


Using  from  Eq.  (5-7)  in  above 


For  a  synchronous  satellite  in  the  ecliptic  plane,  the  major  cau.se  of  deviation 
of  the  satellite  orbit  from  the  ecliptic  plane  is  due  to  rotatic';*  of  the  ascend¬ 
ing  node  of  the  orbit  referred  to  the  equator  [4] . 

In  Fig.  3,  the  normal  to  the  satellite  orbit,  OL,  will  exe-ute  a  conical 
motion  around  the  earth  axis  OZ'.  When  the  orbit  normal  is  OJ. ,  the  orbit 
plane  is  coplanar  with  the  ecliptic  plane.  Tlie  semi-vertical  angle  of  the 
cons  is  q,  which  is  the  inclination  between  the  ecliptic  and  the  equatorial 
plane.  Let  the  ascending  node  of  the  orbit  referred  to  the  equator  and  the 
OY  direction  be  fi.  It  follows  from  Fig.  3,  that  the  inclination  of  the  orbit 
referred  to  the  ecliptic  is 


i  =  2  sin  ^  ^in  q  sin  ^  .  (6-7) 

Similarly,  after  some  manipulations  of  the  rotation  matrix  similar  to  the  treat¬ 
ment  in  Section  IV,  it  can  be  seen  that  the  angle  u  between  OX  and  the  line 
of  nodes  of  the  orbit  with  the  ecliptic  plane  is  given  by 


Nutation  of  satellite  orbit  plane 


_ (l+cos  fl)  cos  n 

(1+cos  cos  ri  +  sin  fi 


(6-10) 


From  Ref.  4,  the  combined  effect  of  the  geopotential  and  lunar  potential 
ia  to  cause  a  precession  of  the  line  of  nodes  between  the  equatorial  plane  and 
the  orbit  plane  at  a  rate  of  dfi  =  —  4.S  deg/year.  The  desired  values  of 
tlie  orbit  elements  referred  to  the  equator  at  the  beginning  of  the  lifetime  are 


i  =  4.65  deg 

SI  =  11.6  deg  . 

With  the  above  initial  elements,  the  elements  at  the  end  of  tiie  five  year 
lifetime  will  be 


(6-11) 


i  =  4.65  deg 

Q  =  —  11.6  deg  . 


(6-12) 


If  we  make  an  allowance  of  0.35  deg  injection  error,  then  the  maximum  inclina¬ 
tion  from  the  ecliptic  at  any  time  will  be  [7]  5.0  deg.  Assuming  S2  =  0  deg 
and  n  =  23.44°  in  Eq.  (6-9)  and  (6-10),  Eq.  (6-6)  gives 


Aii/(t)  =  5a  <0.034  sin 


/  4irt 

135572 


+  2u)|  +  0.008  cos 


/  4TTt 


(6-13) 


The  errci  given  by  Eq.  (6-13)  is  of  a  long  tern  and  can  be  allowed  to  build  up 
and  then  corrected  by  resetting  the  synthesizer.  It  will  be  quite  adequate  to 
reset  the  synthesizer  about  three  or  four  times  over  the  fiw-year  lifetime  of 
the  satellite. 
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VII.  GENERATION  OF  THE  SOLAR  EPHEMERIS 

It  was  seen  in  Section  V  that  the  maximum  error  due  to  discrete  sampling 
of  at  half  day  intervals  is  about  0.0084  deg  and  the  maximum  differ¬ 
ence  between  two  consecutive  is 


Also  from  (3-8)  and  (5-5) 


=  0.017  deg 


(7-1) 


^  1.92  deg 


(7-2) 


It  is  seen  from  above  that  it  will  require  a  much  smaller  memorv'  to  store  the 
difference  Ai;>j(il),  than  to  store  the  function  i/'j^(£.)  itself.  The  function 
ij;j(2.)  can  be  generated  from  A\//(£)  by  using 


<J'j(k)  =  2  Aij/jdc)  +  , 


(7-3) 


where  Ai);(l0  =  'l'(k)  -  >|^(k  -  1)  ,  and  ij^(o)  is  the  function  (|;j^(o)  at  the 
initial  sampling  instant.  The  above  qjeration  is  performed  by  an  adder.  To 
generate  the  ephemeris  for  an  inclined  orbit,  using  Eq.  (5-7),  the  function 
ip,  (k)  has  to  be  read  out  at  proper  phase  and  frequency  and  multiplied  by  the 

factor  ~  .  This  required  multiplication  can  be  performed  by  controlling  the 

“  1 
rate  of  summation.  For  exanple,  multipliers  between  zero  to  j  can  be  obtained 

by  a  jimp  summation 


(7-4) 


where  Int(s)  is  the  integral  part  of  "s"  and  6- -  = 

1  ’  0  l=j 

for  multiplieis  between  y  to  1,  a  skip  simaiation  gives 


Similarly 
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k 

^’i'l  (j) 
j=o 

where 


Above  relations  are  strictly  valid  only  rdien  the  sanpling  interval  is  small 
as  cOTipared  to  the  period  of  the  function  (jf^Ck),  and  u  is  a  small  number. 
IVhen  u  is  a  large  integer  the  mu’  .iplier  is  very  small  and  the  inaccuracies 
in  using  (7-4)  and  (7-5)  for  the  inclinatirai  con'ection  term  are  not  significant 
For  inclinations  ip  to  14.8  deg  only  junp-summation  is  required.  Since  the 
number  u  can  be  only  an  integer,  therefore,  the  inclination  factor  I(i)  can 
be  selected  only  at  discrete  intervals.  The  maximum  error  due  to  discrete  iCi) 
is  a  function  of  the  inclination  i  and  it  is  given  by 


U-sCk) 


(7-6) 


The  multiplier  required  for  various  inclinations  and  the  maxinum  error  due  to 
iuip  summation  in  the  vicinity  of  this  inclination  is  given  in  Table  IV. 

The  error  shown  in  Table  IV  has  a  long  term  period  of  one-half  year.  The 
short  term  error  due  to  the  above  scheme  of  generating  the  inclination  correc¬ 
tion  term  depends  directly  on  the  accuracy  of  the  storing  of  the  differences 

Aii/j(k)  . 
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VIII.  DESIGN  OF  THE  EPHEMERIS  SYNTHESIZER 

The  least  significant  bit  for  the  storage  of  the  differences  is 

chosen  to  be  0.00275  deg.  This  makes  the  synthesizer  compatible  with  the- 
LES-8/9  earth  sensor  which  has  a  least  significant  bit  of  0.011  deg.  With 
the  above  choice  the  error  due  to  the  discrete  representation  of  angle  will  be 
less  than  0.001375  deg.  This  contributes  directly  to  the  short  term  erro"’- 
of  the  synthesizer.  The  nunibei-  representation  is  in  binary  two's  con5)lement. 
The  memory  necessary  for  storage  of  is  a  730  words  of  4  bits, 

including  the  sign  bit.  A  twelve  bit  adder  is  sufficient  to  add  and  store  the 
ephemeris.  The  various  design  features  of  the  ephemeris  synthesizer  are  sunmar 
ized  in  the  Appendix  A.  Appendix  B  illustrates  some  ephemerides  synthesized 
by  a  breadboard  version  of  the  syntliesizer. 

A  block  diagram  of  the  synthesizer  is  shown  in  Fig.  4.  768  word  by 

4  bit  read-only  memory  is  used  for  storing  the  function  (k) .  The  memory' 
power  is  strobed.  l.SW  instantaneous  ii^mory  power  is  stored  in  a  capacitor; 
therefore,  the  maxinura  pov.'er  drained  fran  the  supply  when  reading  the  memory 

is  only  0.3  nW.  The  synthesizer  is  provided  with  the  capability  of  overriding 

the  part  of  the  logic  which  generates  the  ^  ip,  C2k  +  k  )  term  and  generating 

<1  X  o 

only  the  U!j,(k)  term. 
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APPENDIX  A 

DESIGN  SPECIFICATIONS  OF  THE  SOLAR  EPHEMERIS  SYNTHESIZER 


Satellite  orbit 

Circular  synchronous* 

0  to  10  deg  inclination  with  ecliptic 

Satellite  attitude 

The  ephemerides  are  valid  for  a  satellite  having 
no  attitude  errors*. 

Short  term  error 

<  0.001  deg 

Long  term  error 

0.050  deg 

Resettings  (conmand) 

Four 

required 

(1)  Time  of  resetting  12  bit 

(2)  Ascending  node  12  bit 

(3)  Inclination  8  bit 

(4)  Initial  ephemeris  12  bit 

Outputs 

Cl)  12  bit  solar  ephemeris  in  two's  conplement 
(2)  Ephemeris  compute  signal  --  indicates  time 
during  which  the  ephemeris  is  computed 

Other  outputs  available  for  system  check  and 
testing 

Inputs 

(1)  100  kHz  clock 

(2)  10  kHz  clock 

(3)  1/8  day  clock 

Max.  power  required 

700  miV  at  +  5V  ±  10% 

Size 

0.6"  X  4-3/4"  X  6-1/4" 

Weight 

0.72  lbs. 

Additional  corrections  required  for  eccentricity  and  attitude  errors  are 
done  by  using  the  output  of  two  sun  transit  sensors  and  the  infrared  earth 
sensors.  The  computation  of  these  corrections  are  done  in  the  longitude  logic. 
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APPENDIX  B 


EXAMPLES  OF  SOLAR  EPHENCRIS  SYNTHESIZED  FOR  VARIOUS  ORBITS 


In  the  next  few  pages,  some  exanples  of  the  solar  ephemeris  synthesized 
for  various  orbits  are  shown.  The  synthesizer  is  a  mechanization  of  the 
Eq.  (5-7). 


’J'Ct) 


'^'l 


'  Airt 
,365.27 


+  2a) 


“) 


where, 

I 


(B-n 


i  =  inclination  of  the  orbit  with  the  ecliptic  plane 

0)  -  angle  measured  along  the  ecliptic  between  the  direction  of 
perigee  and  the  ascending  node  of  the  orbit. 

t  =  time  measured  frcan  the  instant  of  the  perigee  passage. 


27 


AZIMUTH  (deg) 


Fig.  B-1.  Solar  epheineris  for  a  synchronous  satellite.  No  attitude 
errors,  i  =  0.00°,  o)  =  0°,  e  =  o,  a  =  Rgyj,- 
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AZIMUTH  (deg) 


AZIMUTH  (deg) 


